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Abstract 

Nonperturbative exact solutions are allowed for quantum integrable models in one space-dimension. 
Going beyond this class we propose an alternative Lax matrix approach, exploiting the hidden multi- 
space-time concept in integrable systems and construct a novel nonlinear Schrodinger quantum field 
model in quasi-two dimensions. An intriguing field commutator is discovered, confirming the integra- 
bility of the model and yielding its exact Bethe ansatz solution with rich scattering and bound-state 
properties. The universality of the scheme is expected to cover diverse models, opening up a new 
direction in the field. 
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1 Introduction and Motivation 

A large number of quantum models in one space-dimension (ID) admits exact nonperturbative solutions, 
in spite of their nonlinear interaction. This exclusive class of models, which also includes field models, 
constitute the family of quantum integrable (QI) systems |U [21 El 0, E] with extraordinary properties, 
like association with a quantum Lax and a quantum R matrix, possessing rich underlying algebraic 
structures to satisfy the quantum Yang-Baxter equation (QYBE), existence of a commuting set of 
conserved operators with an exact solution of their eigenvalue problem (EVP), etc. This exact method of 
solution, known as the Bethe ansatz (BA), was pioneered by Bethe way back in 1931 [6j and generalized 
later to algebraic BA mmmm- These QI systems defined in 1 + 1-dimensions, include a wide variety of 
models, e.g. isotropic [6] and anisotropic M quantum spin-^ chains, 5 and 5 -function Bose PHO] and 
anyon pTl 12] gases, nonlinear Schrodinger (NLS) field HI ;13] and lattice [3] model, relativistic m and 
nonrelativistic m Toda chain, t-J |16] and Hubbard [HI [18] model, Gaudin model m, derivative NLS 
[20 ], sine-Gordon [2T] and Liuoville [22] model, etc. The algebraic structures underlying these models are 
also rich and diverse, which include canonical, bosonic, fermionic, anyonic and spin algebras, quantum 
oscillator and quantum group algebras etc., having inherent Hopf algebra properties [23] . However, it is 
important to note, that among this diversity there is a deep unity, revealing that all known QI models, 
we are interested in, are realizable from a single ancestor Lax matrix or its q-deformation [3] El] . At the 
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same time, the diverse algebras underlying these integrable models are also reducible from the ancestor 
algebra or its quantum-deformation [24]. There is a separate class of models with long-range interactions 
[2511261 [27], which although are solvable quantum many body systems, exhibit different properties than 
those listed above and will not be discussed here. The ancestor model scheme, though a significant 
achievement in unifying and generating integrable models, seems to be also an apparent disappointment, 
since it looks like a no-go theorem, allowing no construction of new integrable models beyond the known 
ancestor model. Moreover, since the ancestor model and hence all QI models as its descendants, are 
defined in ID, it apparently excludes any construction of integrable quantum models in higher space 
dimensions. 2D Ivitaev models [28[l29j, belonging to a different class, are possibly the only exception. 

Therefore for a breakthrough, we look for new ideas and observe, that the rational ancestor Lax matrix 
depends on the spectral parameter A only linearly, while its q-deformation depends on its trigonometric 
functions [3j [21. Consequently, all quantum Lax matrices of known integrable models, since realized 
from the ancestor model, depend also linearly (for the rational class), or trigonometrically on A (for 
q-deformed class). For going beyond the prescribed form of the ancestor model, we search for a Lax 
matrix with higher scaling (or length) dimension linked to the integrable hierarchy and for introducing 
extra space dimension, exploit the concept of multi-space-time [H EU x n ,t n , n = 1,2,3,... , hidden 
in integrable systems. For a concrete application we confine to the n = 2 space case and propose an 
alternative Lax matrix approach with A 2 dependence, focusing on the NLS field model as an example. 
It is quite surprising, that though such higher order Lax matrices (with higher order poles) are well 
known in the context of classical integrable systems, they have never been used, as far as we know, in 
the construction of quantum models. Note that, taking x = aq, y = x 2 , t = t;$ in the NLS hierarchy, 
would result to the inclusion of an extra space-dimension y, apart from x and the construction of a novel 
quasi-(2 + 1) dimensional NLS quantum field model, involving a scalar field q(x,y,t ) and its conjugate 
. For confirming the complete integrability of the model, one needs to show the mutual commutativity 
of all its conserved operators, which is guaranteed when the associated Lax matrix satisfies the QYBE. 
However, this task for the present Lax matrix turns out to be the most difficult one, since the commutation 
relations (CR) for the basic fields, known for the existing QI models fail here, due to significantly different 
structure of our Lax matrix and its higher A dependence. Moreover, we can no longer seek the guidance 
of the ancestor algebra [ 23] , since we have gone beyond the known ancestor model. Fortunately, we could 
discover intriguing algebraic relations for our basic quantum fields, which solve the required QYBE with 
the known rational R-matrix. Since the QYBE not only proves the integrability of a quantum model, but 
also gives the CR between the generator of conserved operators and the generalized creation operator, 
we can go ahead with the application of the algebraic BA to our quasi-2D quantum field model and solve 
exactly the EVP for all its conserved operators including the Hamiltonian. Many particle scattering 
and bound states differ considerably from the known result for the ID NLS model. The bound states, 
corresponding to a complex solution for the particle momentum, are found to exhibit unusual properties 
with a variable stability region, dependent on the particle number, coupling constant and the average 
particle momentum. 

2 Quantum integrable models as descendants 

QI models are associated with a discretized quantum Lax matrix A), the operator elements of which, 
for ensuring the integrability of the model, must satisfy certain algebraic relations, which are expressed 
in a compact matrix form through the QYBE 

R( A - fj) U j ( A) ® = U j (n) ® U j (\)R{\ - fj), (1) 

at each lattice site j = 1,2,... IV, together with an ultralocality condition 

[C/ 2 (A)®C/ fc (/i)]=0, j + k. (2) 
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Individual Lax matrices, each representing a particular integrable model, differ substantially in their 
structure, content, nature of the basic fields and underlying algebras, whereas the quantum R-matrix, 
appearing in the QYBE as structure constants, remains the same for all models belonging to the same 
class and therefore can be of only three types: rational, trigonometric and elliptic. However, in spite of 
widely different Lax matrices linked to the rich variety of known QI models, they are in fact realizable 
from a single rational ancestor Lax matrix or its q-deformed trigonometric form [24) . We will not be 
concerned here with the elliptic models, which are anyway few in number. The rational ancestor Lax 
matrix taken in the form 


U r Aric(^) — 


Cl(A + S 3 )+C 2 , s \ 

s + , c 3 (A - s 3 ) - c 4 J ’ 


(3) 


satisfies the QYBE with the well known rational R-matrix [T], due to its underlying generalized spin 
algebra 

[s“,s + ] = 2mV + m“, [s 3 ,s ± ] = ±sj = (4) 

where m + = cic 3 ,m _ = C 2 C 3 + c 4 c 4 , with CjS as Casimir operators or constant parameters admitting 
zero values, and is capable of generating the known quantum integrable models of the rational class. The 
rational quantum R (A — /u) matrix in its 4 x 4 matrix representation may be defined through its nontrivial 
elements as 


R \} = J ?22 = a(A — fj) = A — fj, + ia, 

Ft -21 = Ft -12 = — p) = A — n, i ?22 = R\\ = c = ia, (5) 

while the trigonometric case has q-deformed elements: a = sinh(A— b = sinh(A— /L), c = sinh(m). 
The representative Lax matrices of known QI models of the rational class can be recovered from the 
rational ancestor model ©■ We present below a few of such examples to illuminate the situation. A 
general form for the Lax operators, which can be realized through a bosonic representation from ([3]) was 
proposed earlier [3] (Ch. VIII.4). 

2.1 Generation of rational models 

xxx-spin chain [7]: The Lax matrix may be reduced from the ancestor matrix © at c 4 = c 3 = 
1 , c 2 = c 4 = 0 , giving m + = 1 , rn~ = 0 , which transforms ancestor algebra (|4|) to the spin algebra for 
Pauli matrices. 

Lattice NLS model [3]: The Lax matrix may be obtained from © at the above parameter values, 
by mapping spin operators through the Holstein-Primakov transformation to the bosonic operators: 

Mr Qk\ = 5 Jk t 

NLS field model [lj: The Lax matrix may be recovered from its lattice version at the field limit, giving 
the simple familiar form 

\) < 6 > 

with bosonic field CR: [< 7 ( 2 ;), q'(x')] = 5(x — x'). 

Toda chain [Mj: Lax matrix may be obtained from © at the parameter choice c 4 = 1, c 2 = c 3 = c 4 = 
0 , resulting both = 0, with generators of the reduced algebra realized through canonical variables 
[QjiPk] = djk- 

The rest of the QI models of the rational class, like xxx-Gaudin chain, tJ and Hubbard model 
etc. can also be covered by rational ancestor model ©, employing limiting procedures, higher rank 
representations, fermionic realizations etc., details of which we skip. 
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2.2 Trigonometric models 

Similarly QI models belonging to the trigonometric class, e.g. xxz spin chain, relativistic Toda chain, 
sine-Gordon model, Liuoville model, derivative NLS model etc., are derivable from their representative 
Lax matrices, which in turn can be generated from a single trigonometric ancestor Lax matrix. This 
ancestor matrix is a q-deformation of ([3]) and satisfies the QYBE with the trigonometric i?-matrix, due 
to its underlying generalized quantum group algebra. The details, which we omit here, can be found in 

m- 


3 Novel quasi-2D NLS model 

Since the known quantum Lax matrices (including ([ 6 ]) ), as discussed above, inherit their properties from 
the ancestor models, all of them depend on the spectral parameter A linearly (for rational models) as 
in m, or on sinA, cosA functions (for trigonometric models) and are defined in one-dimensional space. 
A prominent example of rational models is the (1 + l)-dimensional NLS field model associated with the 
Lax matrix ©. Therefore for going beyond the known models and introducing extra dimensions, we 
look into the background concept of multi-space-time dimension {x n }, {t n }, n = 1, 2,... N, hidden in the 
theory of integrable systems [30l [31] . In this formulation of multi-dimension one can define multiple Lax 
equations of the form 


T x„ = U( n -)T, T tn — V( n )T, n — 1,2,... N, (7) 

(here and what follows we denote partial derivatives as subscripts, as a short-hand notation) where 
T = T(X,q ) is the monodromy matrix dependent on the field q = q({x n }, {i n }), defined in multi-space¬ 
time and the generators of the infinitesimal space-time translation Ufa)( A), VLHA) are the space and time 
Lax operators, respectively. However, since the set of Lax equations 0 is a overdetermined system, its 
compatibility conditions (equality of mixed derivatives) would lead to the pairing between any two Lax 
matrices: {Ufa), VLo), (Ufa), Ur m )), (V/ n ), Vr m )), due to symmetry among the variable. Consequently, this 
would lead to the flatness condition among each Lax pair, generating a series of classically integrable 
hierarchal equations as 


dtmUfa) 9 Xn Vfa{) + \Ufa)Vfa x )\ 0, 

dxmUfa) - d Xn U( m) + [Ufa)U( m )\ =0, n,m = 1,2,... N, (8) 

etc. and similarly with other pairs. A possible reduction U n = —V n may be introduced due to the 
simplified condition dt n Ufa) — d Xn Vfa\ = 0. The hierarchal equations ([8]) represent integrable systems in 
(1 + l)-dimensions, 2-dimensions or in quasi-higher dimensions. 

We intend to use this concept of the hierarchy of multiple space-times, embedded in integrable systems, 
for constructing quantum integrable models in quasi-(2+l) dimensions, restricting to the case q(x\,X2,t^), 
involving two space and one time Lax operators (Ufa), Ufa), V( 3 )), linked to space x = x\,y = X 2 and time 
t = t 3 variables. 


3.1 Alternative Lax matrix 

For a concrete application, we consider the NLS family of field models, which belongs to the rational 
class and choose our quantum Lax matrix as its next hierarchy: 


U ( 2 ) (A) 


2A 2 - q^q 2A q - iq x A 

-2A- iq\ —2A 2 + q^q ) 


(9) 
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It is interesting to compare the structure of Lax matrix © with that of the well known NLS model ([ 6 |), 
to note the crucial differences, that the matrix elements of ([9]) depend on the spectral parameter up to 
A 2 (double pole) and involve field operators q, q ', q x , q' xl defined in quasi -(2 + 1 ) dimensions: (x,y,t). It 
needs to be mentioned, that such higher order Lax matrices like © appearing in the integrable hierarchy 
are usually taken to be independent entries, not constructed solely out of lower order Lax matrices, 
though there are formalisms to connect them in an involved way using classical r-matrix [32]. However 
no quantum extension of this method is available and it is also not clear whether the general L-operators 
of [3, 2i] can be used for this purpose. Therefore, leaving aside the question about the possibility of 
constructing ([9]) from more fundamental Lax operators, we start directly with this higher order Lax 
operator for constructing our quantum model. 

We emphasize again that, in spite of the familiarity of Lax matrix ([9]) in classical integrable hierarchy, 
such higher-pole Lax operators have been ignored so far in the context of quantum integrable models. 
Note that, in the hierarchal equations we consider the pairing (Ur 2 -\,V^) with x = x±,y = x 2 ,t = t%, 
for constructing quasi -(2 + 1 ) dimensional model, which can also be reached by a combination of the 
Lax pairs (Z7m, H( 2 )) in 2-dimensions with x = x\ , y = x 2 and (U^,V^) in (1 + l)-dimensions with 
x = xi,t = t^. Note that, the x-shift space Lax operator here can be given by Um in ©, associated 
with the space Lax operator of the standard NLS model, while the y-shift space Lax operator is given by 
17(2) i n © and the f-shift time Lax operator by V( 3 ), representing a higher order A 3 dependence (cubic 
pole) form, which we do not reproduce here. However, for constructing our quantum model and exactly 
solving it through algebraic Bethe ansatz we would need only the quantum Lax operator U( 2 \. 

We introduce here the notion of scaling or length dimension, which is an useful concept in analyzing 
higher order Lax operators in multi-space-time dimensions. One defines a scaling dimension [L~ l ] = 1 
for length L. Therefore from © we get [U^] = n, since [d Xn ] = n (similarly for VL) ) and consequently 
[[/(!)] = 1, since each term in © has scaling dimension [A] = [g] = 1. Similarly, [17(2)] = 2, since in © 
[A 2 ] = [|g| 2 ] = [q x \ = 2 , and [V( 3 )] = 3, etc. 

3.2 Quantum integrability through Yang-Baxter equation 

In dealing with quantum field models one has to lattice regularize the Lax operators first to avoid short- 
distance singularities [T]. Therefore, our intention is to show, that the discretized Lax matrix along the 
y-direction: W = I + AI7( 2 )(A, qj). where qj = q(x,y = j,t ), with lattice constant A —> 0, does satisfy 
the QYBE © with the rational .R-matrix ©. However, this becomes a highly involved problem, since 
due to more complicated structure of the present Lax matrix ©, ten out of total 16 relations of the 
4x4 matrix QYBE remain nontrivial, all of which are to be satisfied with a suitable field CR. Compare 
this situation with the known ID NLS case 121, where due to much simpler form of the Lax matrix 
?7(i) ©, only two nontrivial relations in the QYBE survive, which can be solved successfully using the 
bosonic field CR. However, we realize that, no algebraic relations, including the bosonic CR, appearing 
in the existing integrable models would work here, since the choice of C/( 2 ) has taken us beyond the scope 
of the known ancestor models and the associated algebras. Moreover, the CRs for the field now have to 
be sought for along the extra direction y , that has been included in the system. Therefore we look for 
some innovative algebraic relations for the basic quantum fields to be consistent with the QYBE, linked 
to the present Lax matrix ©. Fortunately, we find a new set of such relations for our quasi-2D fields as 


[q{x,y,t), qi{x,y',t)] = -2ia 6(y - y'), [ q x (x,y,t ), q\x, y, Q] = 2iaS(y - y'), (10) 

[q{x,y,t), q ] (x,y',t)} = 0, (11) 

(or their discretized version (1231) 1 together with their hermitian conjugates. Note that CRs (HOI) (or 
(1231) 1. exhibiting an asymmetry in space variables are fundamentally new relations, different from known 
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relations like canonical, bosonic etc. It may be observed that, the form of CR (1101) may be linked to 
the quadratic space-velocity term : q\q x : appearing in Hamiltonian (fT3j) . Application of these algebraic 
relations (1 10 1 11 1) satisfies miraculously all ten nontrivial equations appearing in the QYBE, involving the 
discretized Lax matrix U 3 (A) , up to order 0(A) (see 1)2411251 ) in Appendix for details). This is however 
enough for proving the integrability of the field models, obtained at the limit A —>• 0. It is remarkable, 
that in spite of the presence of a x-derivative term, new CRs (jlOD satisfy the necessary ultralocality 
condition ([21). This is because not x but y is the relevant space direction here, where the fields commute 
at space separated points along j —> y, reflecting the quasi-2D nature of our model. 

Therefore, since due to mm (or its discretized version (123]) ) , the lattice regularized quantum 
Lax operator U 3 ( A) constructed from ([9]) satisfies the QYBE (jX|) for the rational i?-matrix, together 
with the ultralocality condition ([2|), the transition matrix for our model, defined for IV-lattice sites: 
T( A) = nf=i U 3 (\), must also satisfy the QYBE [TJ 

R(X - fj) T(A) ® T(y) = T(/x) <g> T(X)R(X - fj), T( A) = ^ ) > ( 12 ) 

with the same R(X — / u)-matrix. This happens due to the coproduct property of the underlying Hopf 
algebra [23], which keeps an algebra invariant under its tensor product. This global QYBE (|12l) serves 
two important purposes. First, it proves the quantum integrability of the model by showing the mutual 
commutativity of all conserved operators. Second, it derives the commutation relations between the 
operator elements of T(A), which are used for the exact algebraic Bethe ansatz solution of the EVP. 

In more details: multiplying QYBE (1121) from left by i? _1 , taking the trace from both sides and using 
the property of cyclic rotation of matrices under the trace, one can show that r(A) = trace T(A) commutes: 
[t(A), t(p)] = 0. This in turn leads to the Liuoville integrability condition: [Cj, Ck ] = 0, j, A: = 1,2,-- ■, 
since the conserved set of local operators are generated from lnr(A) = Yhj CjX~ 3 , through expansion in 
the spectral parameter A. Following this construction and exploiting the explicit form of the Lax matrix 
([9]), we can derive, in principle, all conserved operators Cj, j = 1,2 ,... for our model. Skipping the 
details, which can be found for the classical case in m, we present here only the x-shift Hamiltonian as 
H (x) = C 2 : 

H {x) = J dy : {iq ] q y + qUx + <? t V) : (13) 

and the t-shift Hamiltonian as H = C 4 : 

H = J dy : (: iq\q xy + q\q y + i{q ] q){q ] q y ~ Qyq) 

~ 2 {q^q) q\q x + q ]2 q 2 x + qfq 2 ) :, (14) 

which we take as our model Hamiltonian, we are interested in. Notice the quasi (2 + 1) dimensional 
nature of Hamiltonian ([T4]) . since though it involves both x and y derivatives of the field : q x (x,y,t ) and 
q y (x, y, t), the volume integral is taken only along y. Asymmetry in the appearance of space derivatives is 
also explicit. However, at the same time an operator with double space volume integral: H = -^ f_ L dx H 
in a strip of x £ [-L, L] is also conserved in time, since dtH = 0. 

4 Algebraic Bethe ansatz for the eigenvalue problem 

Since T(A) satisfies the QYBE ( 1121 ) with the rational R- matrix, we can follow the procedure for the 
algebraic BA, close to the formulation for the ID quantum NLS model mm • As we have discussed 
above, r(A) = traceT(A) = A(A) + A (A) is linked to the generator of the conserved operators Cj, j = 
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1,2,..., including the Hamiltonian (1141) . The off-diagonal elements of T(A) (|12D : B( A) and L>f(A), on 
the other hand, can be considered as generalized creation and annihilation operators, respectively. For 
solving the EVP for all conserved operators: Cj \ M >= Cj hf> \M >, j = 1,2,... simultaneously, we 
construct exact M-particle Bethe state |M >= B(yi)B(y 2 ) • • • >, on a pseudo-vacuum |0 > 

with the property B\yj) |0 >= 0, A(A)|0 >= ao(A)|0 >, and aim to solve the EVP: r(A)|M >= 
Am(A,/xi,/X 2 , • • • ,Hm)\M >, with exact eigenvalues In A M (A, {//„}) = Ej cf I \{y a })X~ j . 


4.1 Exact solution for quasi 2D quantum field model 


For obtaining the final result for our quantum NLS field model on infinite space interval we follow the 
formulation of [I] for the ID NLS field model, though adopted here for the transverse dimension y and 
higher order Lax operator ([9]). We switch over to the field limit: A — > 0 with total lattice site N — > oo and 
then take the interval L = N A —> oo, assuming vanishing of the field qj —> 0, at j —> oo, compatible with 
the natural boundary condition of having the vacuum state at space infinities, yielding the asymptotic 
Lax matrix W (A)|j_>.oo = Ho (A) = 2iX 2 a 3 . Therefore, we have to shift over to the field transition matrix 
defined as 

T f ( A) = U^ N T{ A) Uq N , N —► oo, (15) 

and for further construction introduce V(X,y) = Fo(A) (g) Uo(y), W(X,y) = (tV(A) ® W (y))j^oo. We 
may check from the QYBE dTJ that W satisfies the relation R (A — y)W (A, y) = W(y, \)R(X — y), using 
which we can derive from QYBE (1121) . that the field transition matrix (|15D also satisfies the QYBE 

i?o(A, y) T f {\) (g> Tf(y) = T f (y ) <g> T f (\)R 0 (\, y), (16) 

but with a transformed i?-matrix: 

R 0 = S(y, X)R(X - y)S( A, y), S( A, y) = W~ N { A, y)V N { A, y), N ^ oo, (17) 

where R(X — y) is the original rational R —matrix ([5]) (see [lj for similar details on ID NLS model). 

Based on the above formulation and using the field operator products: <Rq 3 ^ x = —2 ij^, <i ] X qj = 0, 
at j —> oo, compatible with the field CR, we can calculate explicitly the relevant objects needed for our 
field model. In particular, the central 2x2 block W c for matrix W turns out to be 

^(A,m) = / + AM(A, (1 )( (A __ (A 0 _„))- (IB) 

with an intriguing factorization of its spectral dependence by a prefactor M(X,y ) = 2(A + y), which is 
the key reason behind the success of the exact algebraic Bethe ansatz solution for our field model, in spite 
of the more complicated form of its Lax operator. Note, that since our model shares the same rational 
i?-matrix with the known NLS case (though having different Lax operators), the present result coincides 
in part with that of the ID NLS model mug, though only formally. On the other hand, the transformed 
Rq matrix, relevant for the field model, depends on the corresponding asymptotic Lax matrix and its 
product through matrix S(X,y). Therefore, since Lax matrix @ for our model is more complicated, 
compared to d 6 j) for the ID NLS model, our final result shows intriguing differences from the known NLS 
result, which we highlight below. 

For constructing Rq using definition (1171) . we have to construction first matrix S(X,y), taking proper 
limit of W~ N at L —> oo using (|18l) . Through some algebraic manipulations, which are skipped here, 
we finally arrive at the field limit to a simple expression for R$ matrix, expressed through its nontrivial 
elements as 


Ru = = a(A - y), R l 2 ( = b(X - y), R& = Rff = 0, 

,,2 




cr 


+ 


a 2 7r 


X — y M(X,y) 


S(X-y), 


(19) 
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where M(A, n) = 2(A + /x), the terms a{ A — /x), 6 (A — /x) are as in ([5]) and the 6 (A — //) term vanishes 
at A / /*. It is interesting to compare i?o-nratrix (I19p . relevant for the field models, with the original 
R-matrix (j5|) . Now from QYBE (1161) linked to field models, we can derive for our model the required 
CR between the operator elements of Tf, using Rq matrix (1191) . In particular, we get the commutation 
relation 

A - fij) + M °^ ^ ~ A), (20) 

where fj = - E. — . Note that the prefactor fij) = appearing in the above CR bears 

the imprint of the A 2 dependence of our Lax matrix and is absent in such relations in the standard NLS 
model. At A ^ jij however, when the singular term vanishes, the relation coincides formally with the 
known NLS case. Using this result and the property of the vacuum state |0 >, we obtain the exact EVP 
for 

M 

A(X)\M >= Fm\M >, as F M = Wf j {\-ii j ) 

3 

and hence for r(A), which gives finally the exact eigenvalues for all conserved operators. For our model 
Hamiltonian H = C 4 , we obtain the exact energy spectrum Em = Yj 1 Rj? f° r the M-particle scattering 
state, which clearly differs from that of the known NLS model M, though bearing formal similarity 
with the next NLS hierarchy. 

It is remarkable, that in spite of the highly nonlinear field interactions present in the Hamiltonian (1140 . 
the scattering spectrum shows no coupling between individual quasi-particles, mimicking a free-particle 
like scenario. On the other hand, the bound-state or the quantum soliton state, which is obtained for the 
complex string solution for the particle momentum: = fiQ + x^((M + 1 ) — j), where /x 0 is the average 

particle momentum and a is the coupling constant, induces mutual interaction between the particles. 
The corresponding energy spectrum may be given by 

M 

E { M Und) = #e[f>$ s) ) 4 ] = M+ E b (v 0 , a, M), (21) 

3 

where E b is the binding energy of the M > 1-particle bound-state. Recall, that a bound-state becomes 
stable, when its energy is lower than the sum of the individual free-particle energies, which in turn is 
ensured by the negative values of the binding energy: E b . More negative binding energy indicates more 
stable bound-states. For the known NLS model the binding energy mH3] Ef a (a,M) = —gM(Af 2 — 1 ) 
is independent of no and strictly negative, which makes the corresponding bound-states always stable 
with the stability increasing as the particle number M and the coupling constant a increase. 

However, for the present quasi-2D NLS model, the picture differs significantly, producing a fascinating 
bound-state spectrum with intricate stability region. Note, that in the present case the binding energy 

2 2 

E b (no, a, M) = E+ + E-,E- = -^JH[(M 2 - 1) 

E+ = ^M[(^M 4 + 1 - 2 -M(M + ^))], (22) 

contains both negative and positive terms (due to the simple mathematical fact, that in the expression 
EfOxW ) 4 the real term with i 2 = —1 gives negative, while i 4 = +1 strictly positive contribution). 
Therefore, binding energy (122j) may take negative as well as positive values, depending on the parameters 
a, IM) and M, defining a variable domain for the existence of stable bound states (see Fig.l). Note that 
the term E~ = proportional to that of the known NLS model, stabilizes the bound-state, while 

the counter term E + has a destabilizing effect. The graphs show clearly that the stability region of 






Binding Energy 




Figure 1: Binding energy for the quasi-2D NLS model. Figure shows E & (I22j) with increasing particle 
number M for different values of the coupling constant a =0.1 (red), 0.5 (pink), 0.7 (black), 1.0 (blue), 
2.0 (green), with parameter //q fixed at (a) 1.0 and (b) 2.4. 


the bound-state for our model shrinks with the increase of the particle number M as well as with the 
coupling constant a, which is rather anti-intuitive, since for the known NLS model,as seen from E£ ls , 
the bound-state stability always increases with increasing M and a. On the other hand, one can enlarge 
the stability domain in our model by increasing the particle momentum /io, as evident by comparing 
the figures (a) and (b). This feature however can not be matched with the known result of the ID NLS 
model, since its binding-energy is independent of [iq. This shows the intricate nature of the bound-state 
configuration for our model in comparison with known result of the ID NLS model. 

5 Conclusion and outlook 

Going beyond the known form of the existing integrable quantum models in ID, we propose an alter¬ 
native higher order Lax matrix approach, exploiting the concept of multi-space-time dimension hidden 
in integrable systems, and apply it for constructing and solving a novel quasi-2D quantum NLS field 
model. The key to our success in proving the crucial quantum Yang-Baxter equation, which guarantees 
the quantum integrability of the model, is the discovery of a new type of field operator algebra, not 
covered by the existing rules. 

The known ID quantum integrable models satisfying QYBE with rational R-matrix may be realized 
from an ancestor Lax matrix associated with a spin-like algebra, reducible to conventional spin, bosonic, 
or the canonical algebra, related to the existing models. A bosonic realization of this matrix was proposed 
as a general L-operator in Ch. VIII.4 of [3], followed by a theorem, stating that the same L-operator can 
construct a monodromy matrix with arbitrary rational function. However this theorem, claiming only a 
sufficient but not a necessary condition, does not rule out the possibility of alternative L-operators, an 
example of which is provided by the present Lax matrix. 

We stress that, the quantum Lax operator ([9]) with higher scaling dimension and the associated novel 
commutation relations (jlOllllft for the fields in our model are fundamentally different from those used in 
nonrelativistic integrable quantum systems. Unlike the known L-operators belonging to the rational class, 


9 











(}9j) with A 2 spectral dependence and having x-derivative of the field, can not be realized straightforwardly 
from any linear combination of the general ultralocal L-operators proposed earlier [3 El]. Similarly, the 
crucial algebraic structure (HOD (with x-derivative term) is different from known ultralocal algebras and 
evidently can not be generated by combining them. However, a possible construction of such higher- 
order Lax operators as a nonlinear combination (like product) of lower order Lax operators of [3 (21] and 
obtaining the underlying novel algebras from the known ones could be taken up as a challenging future 
problem. 

The dimensionality of the present model with its field q(x. y, t ) needs special focus. In one hand, the 
system shares effectively one-dimensional properties, since it is linked to the ID NLS hierarchy. This 
is also reflected in the energy spectrum of the present model, which is similar to that of the higher 
Hamiltonian in the known NLS hierarchy. On the other hand, the model Hamiltonian (I14D contains 
derivatives of the field q x , q y in both x and y variables and similarly both these variables are involved in 
commutators (I10II11D as well as in the present Lax matrix, defining the model in quasi two-dimensional 
form. Moreover, these 2D structures can not be reduced to ID by ignoring the dependence on the other 
variable. 

Due to quantum integrability of our quasi-2D NLS model, the eigenvalue problem can be solved 
exactly for the commuting set of all its conserved operators, with intriguing result for the many particle 
scattering and bound states. 

It is worth adding that, recently we have constructed a novel quasi 2D quantum Landau-Lifshits model 
belonging also to the rational class (to be reported elsewhere). It is reasonable to assume therefore, that 
such quasi 2D quantum models generated by higher-order Lax operators, are not limited only to the 
present NLS case, but constitute a novel family of quantum integrable systems within the rational class. 
The present approach, general enough for applying to other quasi higher dimensional quantum models, 
could open up a new direction in the theory of quantum integrable systems. It is a challenge to find a 
possible q-deformation of the algebra found here, which could lead to a novel class of quantum algebra, 
while an exact lattice version of the present Lax matrix could unravel a higher-order ancestor Lax operator 
for generating a new family of integrable quantum models. 

6 Appendix 

In QYBE (jT]) with R -matrix ([5]) and discretized version U 3 of the quantum Lax matrix ([9]), out of total 

16 matrix operator relations, except 4 diagonal and 2 extreme off-diagonal terms, all other 10 relations 

Q^j stand nontrivial and their validity needs to be proved using the CR, discretized from (11011 11 j) : 

+ Ot + 

[Qj,Qj, x \ = = ° ( 23 ) 

and their conjugates. 

6.1 QYBE relation for matrix elements 

Using expressions for a(A — y), b(X — y), c of ([5]) and CR (f24l) we may check the validity of 

Qll = a W U (X)U j M - b - c W n (n)U j 12 W) 

= iA(A - n) g(-A[g], qjx ] + 2c) + 0(A 2 ) = 0, (24) 

upto order 0(A 2 ). Similarly one proves the conjugate relations Ql\, Qf{, Q\i and similar relations 

q22 q22 q 12 f )21 

Vl2)L;21> L?22 iL;22- 
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The validity of the remaining two relations can also be proved with the use of the same CR (|23l) : 


Qli = b [U j 12 (A), UP iGu)] + c (^ 22 (A)^'h(m) - U>n(\)U j 22 (p) 

= 2iA 2 (A - fJ.)(/J,[qj,x,qj\ + P<lj, x ,Qj]) + 4iAc(^ 2 - A 2 ) = 0, (25) 

which is valid exactly in all orders of A and similarly for the conjugate relation Q\\. This proves thus 
the validity of all QYBE relations for our quantum quasi-2D NLS field model, associated with the higher 
Lax operator Q and algebraic relations mm, obtained at the limit A —>• 0 . 
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